Exactly solvable Z^Ar-type quantum spin models with long-range interaction 
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We derive the spectra of the Djv-type Calogero (rational) su(m) spin model, including the de- 
generacy factors of all energy levels. By taking the strong coupling limit of this model, in which 
its spin and dynamical degrees of freedom decouple, we compute the exact partition function of the 
su(m) Polychronakos-Frahm spin chain of Djv type. With the help of this partition function we 
study several statistical properties of the chain's spectrum, such as the density of energy levels and 
the distribution of spacings between consecutive levels. 
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I. INTRODUCTION 

Recent studies of quantum integrable dynamical mod- 
els and spin chains with long-range interactions [l|, 0, 
0; 0j S S 0i Si have not only enriched our understand- 
ing of strongly correlated many-particle systems in one 
dimension, but also influenced several branches of math- 
ematics in a significant way. In particular, it is found 
that this class of quantum integrable systems have close 
connections with apparently diverse subjects like gener- 
alized exclusion statistics j3, 0, [l^l, quantum Hall ef- 
fect 11 ill, quantum electric transport in mesoscopic sys- 
tems [l2, fl3{ , random matrix theory , multivariate or- 
thogonal polynomials [isl . [TgI [TtI and Yangian quantum 
groups [13, [H, ■ The interest in quantum integrable 
models with long-range interaction was initiated by a 
seminal work of Calogero [1] , where the exact spectrum of 
an iV-particle system on a line with two-body interactions 
inversely proportional to the square of their distances and 
subject to a confining harmonic potential was computed 
in closed form. An exactly solvable trigonometric variant 
of the rational model introduced by Calogero was pro- 
posed shortly afterwards by Sutherland 2, 3J. The parti- 
cles in this so-called Sutherland model move on a circle, 
with two-body interactions proportional to the inverse 
square of their chord distances. Subsequently, Olshanet- 
sky and Perelomov established the existence of an un- 
derlying root system structure for both the Calogero 
and Sutherland models, and constructed generalizations 
thereof associated with other classical (extended) root 
systems like B^, Cm and BCn 4]. 

In a parallel development, Haldane and Shastry found 
an exactly solvable quantum spin-i chain with long- 
range interactions, whose ground state coincides with the 
U OO limit of Gutzwiller's variational wave function 
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for the Hubbard model, and provides a one-dimensional 
realization of the resonating valence bond state d, 0]. 
The lattice sites of this su(2) Haldane-Shastry (HS) spin 
chain are equally spaced on a circle, all spins interact- 
ing with each other through pairwise exchange interac- 
tions inversely proportional to the square of their chord 
distances. A close relation between the HS chain and 
the su(m) spin generalization of the original (AN-type) 
Sutherland model [U, [S^, |2^, which leads to many 
quantitative predictions, was subsequently established 
through the so-called "freezing trick" 0, US] ■ More pre- 
cisely, it is found that in the strong coupling limit the 
particles in the spin Sutherland model "freeze" at the 
coordinates of the equilibrium position of the scalar part 
of the potential, and the dynamical and spin degrees 
of freedom decouple. The equilibrium coordinates co- 
incide with the equally spaced lattice points of the HS 
spin chain, so that the decoupled spin degrees of freedom 
are governed by the Hamiltonian of the su(to) HS model. 
Moreover, in this freezing limit the conserved quantities 
of the spin Sutherland model immediately yield those of 
the HS spin chain, thereby explaining its complete in- 
tegrability. By applying the freezing trick to the A^- 
type rational Calogero model with spin degrees of free- 
dom, a new integrable spin chain with long-range inter- 
action was constructed in Ref. 0] ■ The sites of this chain 
— commonly known in the literature as the Polychron- 
akos or Polychronakos-Frahm (PF) spin chain — are un- 
equally spaced on a line, and in fact coincide with the 
zeros of the A^-th Hcrmitc polynomial [25j . 

The powerful technique of the freezing trick was sub- 
sequently used to compute the exact partition functions 
of both the su(m) PF spin chain [2^ and the su(m) HS 
chain [131 , the BCn counterparts of these chains [H, |2^ , 
and their supersymmetric extensions [sol. Isil. Is^ . The ex- 
act computation of the partition functions of these quan- 
tum integrable spin chains has opened up the exciting 
possibility of studying various statistical properties of 
their energy spectra. Indeed, it is found that for a large 
number of lattice sites the energy level density of such 
chains follows the Gaussian distribution with a high de- 
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gree of accuracy [27|, [281, [291, l31|, [33, l33| . It has also been 
observed that the distribution of the (normaUzed) spac- 
ings between consecutive energy levels of these chains is 
not of Poisson type, as might be expected in view of a 
well-known conjecture of Berry and Tabor [s^l- An ana- 
lytical expression, which explains the unexpected distri- 
bution of spacings between consecutive energy levels in 
the above mentioned chains, has recently been derived 
using only a few simple properties of their spectra [28j . 

Our aim in this article is first of all to analyze the spec- 
tra of the su(m) spin Calogero model of 13 at type. We 
shall then apply the freezing trick to compute the exact 
partition function of the Dj^ version of the PF spin chain, 
and use this partition function to study various statistical 
properties of the chain's spectrum. It should be stressed 
that, although the Hamiltonian of the D^v-type su(m) 
spin Calogero model can be obtained by setting to zero 
one of the coupling constants of their BCm counterparts, 
this fact does not allow one to find out all physically rel- 
evant properties of the Dm model as a limiting case of 
its BCn version. For example, as will be explained in 
Section [TTl the configuration space of the ZJ^v-type spin 
Calogero model differs quite significantly from its BCm 
counterpart. A more drastic change occurs in the Hilbcrt 
space of the Dm model, which gets "doubled" in com- 
parison with the BCm one. More precisely, the Hilbert 
space of the Dm spin Calogero model can be expressed 
as a direct sum of the Hilbert spaces associated to two 
different BCm models with opposite "chiralities" . These 
remarkable properties of the Dm model indicate that it 
is a "singular limit" of the corresponding BCm model, 
worthy of consideration in its own right. 

The paper is organized as follows. In Section |TT] we in- 
troduce the su(m) spin Calogero model oi Dm type and 
construct its associated (antiferromagnetic) spin chain by 
means of the freezing trick, discussing their relation with 
their BCm counterparts. Section [ni] is devoted to the 
evaluation of the spectrum of the spin Calogero model 
oi Dm type, which is then used to compute in closed 
form the partition function of its associated spin chain 
applying the freezing trick. We also show how to ex- 
press this partition function in terms of those of the PF 
chains of type A and B. In Section IIVI we make use of 
the closed-form expressions for the partition function of 
the PF chain oi Dm type to analyze several statistical 
properties of its spectrum. We show that — as is the case 
with other chains of HS type — the level density follows 
with great accuracy the Gaussian law when the number 
of lattice sites is sufficiently large. We also prove that the 
cumulative distribution of spacings between consecutive 
levels follows the same "square root of a logarithm" law 
obeyed by the PF chain of types A and B and by the 
original HS chain. This provides further confirmation of 
the fact that spin chains of HS type are exceptional inte- 
grable systems from the point of view of the Berry-Tabor 
conjecture. Finally, in Section |V] we outline the general- 
ization of the above results to the ferromagnetic chain 
and its associated spin dynamical model. 



II. THE MODEL 

Since the su(rn) spin Calogero model oi Dm type is 
closely related to its BCm counterpart, we shall start by 
briefly reviewing the latter model, whose Hamiltonian is 
given by [1^ 
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Here the sums run from 1 to (as always hereafter, 
unless otherwise stated), a > 1/2, 6 > 0, e = ±1, 
xfj = Xi zL Xj, = J2i^i^ ^ij is the operator which 
permutes the i-th and j-th spins. Si is the operator re- 
versing the z-th spin, and Sij = SiSjSij. Note that the 
spin operators Sij and Si can be expressed in terms of 
the fundamental su(m) spin generators J" at the site k 
(with the normalization tr{J^J2) — ^S"'') as 



Sij — 

m 



E 

Q = l 



Ji J J , 
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The configuration space of the Hamiltonian ([T]) can be 
taken as one of the Weyl chambers of the BCm root sys- 
tem, i.e., one of the maximal open subsets of on which 
the linear functionals Xi ± Xj and Xi have constant signs. 
We shall choose as configuration space the principal Weyl 
chamber 



C(B) ={x= (xi. 



,XM) 



< Xi < X2 < ■ ■ ■ < Xm} ■ 



(2) 

The spectrum of the BCm spin Calogero model, includ- 
ing the degeneracy factors of all energy levels, has been 
determined by constructing a (non-orthonormal) basis of 
the Hilbcrt space in which the Hamiltonian ^ is trian- 
gular [28] . By setting h = (3a and taking the limit a — > oo 
in the Hamiltonian ([1]) , one can obtain the su(m) PF spin 
chain of BCm type, with Hamiltonian given by 



1 + Sr 



l + S^ 



'E 



H(^) = E 

(3) 

Here /3 is a positive real parameter, and the lattice sites 
S.i can be expressed in terms of the zeros yi of the La- 
guerre polynomial as yi — £,f /2. The exact parti- 

tion function of the spin model ([3]) has also been recently 
computed with the help of freezing trick [28j . 

The Hamiltonian of the su(m) spin Calogero model 
oi Dm type is obtained by setting 6 = in its BCm 
counterpart ([l]), namely 



^--EC + T^ 



Si' 



Sij 



i^j? 



(4) 



3 



As configuration space of the Hamiltonian ([4]) we can take 
again one of the Weyl chambers of the Djv root system. 
For instance, the choice xi < ■ ■ ■ < xn determines all the 
differences Xi — Xj . If we also require that X1+X2 > the 
sign of all the sums Xi + Xj is determined as well. Indeed, 
|a;i| < X2 implies that |a;i| < Xj for all j = 2, . . . , N, so 
that xi + Xj > for j > 1, while the sums Xi + Xj 
with i,j > 2 and i ^ j are clearly positive on account 
of the positivity of Xk with fc > 1. Thus we can take as 
configuration space of H the open set 

C = {x = (xi, . . .,xn) ■■ \xi\ < X2 < ■ ■ ■ < Xn} , (5) 

which is just the principal Weyl chamber of the Dn root 
system. It is interesting to observe that this configuration 
space contains its BCn counterpart ((2|) as a subset. 

The Hamiltonian of the su(r7i) PF chain oi Dm type 
can be obtained from the spin Hamiltonian ([4]) in the 
limit a — > 00 by means of the freezing trick. More pre- 
cisely, since 



H 



dl +a^U + Oia). 



with 



i4r 



4 ' 



(6) 



when the coupling constant a tends to infinity the parti- 
cles in the spin dynamical model ^ concentrate at the 
coordinates of the minimum ^ of the potential U in C. 
From the identity 



a7i(x) . 



(7) 



where 



Hsc = -J2dl + ^r' + aia-l)Y: 



and 



(8) 



1 + 



1 + s,, 

i^Xi -|- Xj ^'^ 



00 the internal degrees 



it follows that in the limit a 
of freedom of H are governed by the Hamiltonian H — 
H(^), explicitly given by 



1 + s.. 



1 + s., 



(9) 



Equation ^ is the Hamiltonian of the (antiferromag- 
netic) su(to) PF chain of type, whose sites £,i are 
the coordinates of the unique minimum ^ of the scalar 
potential ^ in the open set (O. The existence of this 
minimum follows from the fact that U tends to -f 00 on 
the boundary of C and as r — > 00, and its uniqueness 



was established in Ref. '35'| by expressing the potential 
U in terms of the logarithm of the ground state p of the 
scalar Dn Calogero model Hsc, given by 



p(x) = e 4 ' 



ni 

i<j 



(10) 



As shown in the latter reference, the sites coincide with 
the coordinates of the (unique) critical point of log p in 
C, and therefore satisfy the nonlinear system 



ME 



1 



0. 



The numbers cannot be all different from zero, since 
in that case we would obtain the contradiction 



o = E 



1 



■2 

'J 



N 

T 



N 

T 



Hence = for some i, and since (^1, 
we must have 

ei = o, 



, ^n) lies in C 



(11a) 



while the remaining iV — 1 sites should satisfy the condi- 
tion 



E 



1 1 



£2 _ f2 4 



z = 2,...,iV. 



(lib) 



Substituting Eq. (|lla|) into (jllbp one obtains 



N 

E 72 



(12) 



It is interesting to compare the above condition satis- 
fied by the nonzero ^i's with the relation 



2y. 



(13) 



obeyed by the zeros yi of the Laguerre polynomial 



^ It is evident that Eq. ^ reduces to Eq. (fTS]) 
when M = iV - 1, /? = 2 and = £,f^i/2. We there- 
fore conclude that the sites £,2 < ■ ■ ■ < £.n are expressed 
in terms of the iV — 1 zeros yi < ■ ■ ■ < yn-i of the 
Laguerre polynomial L\j_-y by = \J2yi-\. On the 
other hand, it has already been mentioned that the lat- 
tice sites of the PF model of BCn type ([3]) are expressed 
in terms of the zeros of the Laguerre polynomial L^"^ 
by yi — Since the potential U in Eq. ^ is ob- 

tained from its BCn counterpart in the limit /? ^ 0, we 
could also have argued that the lattice site £i of the Dn- 
type PF model is the square root of twice the i-th zero 
of L'^ for % — 1, . . . ,N. The equivalence of both char- 
acterizations is substantiated by the well-known identity 



L],\y) = -yL],_,iy)/N, cf. i35,]. 
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It is worth pointing out that, even though the lattice 
sites of the BC^v-type PF chain coincide with their ZJ^v 
counterparts in the Hmit /3 ^ 0, the Hamiltonian ([3]) of 
the PF chain of BCn type does not reduce to its Dn 
variant ([9]) in the same hmit. To estabhsh this fact, note 
first that all roots of the equation L^^{y) ~ except 
the smallest one tend to a finite nonzero value in the 
limit /? ^ 0. As a result, terms like /3(1 — eSi)/^f, which 
appear in the r.h.s. of Eq. ([3]), vanish for i = 2, . . . , A^. 
We next examine the behavior of the smallest root f i of 
the equation L^^{y) — 0. It can be shown [2^ that the 
zeros of satisfies the relation 



Ref. [231 we introduce the auxiliary operator 



^ 1 



1 

yi 



(14) 



Since the l.h.s. of this equation vanishes in the limit /3 
0, the r.h.s. yields lim/3_^o(2/?/Ci) = N. Substituting this 
limiting value in Eq. ([3]) we find that the Hamiltonians 
of the BCn- and Z^Ar-type PF spin chains are related by 



hm 



(15) 



It is interesting to observe that the second term in the 
r.h.s. of the previous equation may be interpreted as an 
"impurity" interaction at the left end of the BCn spin 
chain. 



III. SPECTRUM AND PARTITION FUNCTION 

We shall start by deriving the spectra and partition 
functions of the DAr-type su(r7i) spin Calogero model (|4]) 
and its scalar counterpart ((8]). Since the spin and dynam- 
ical degrees of freedom of the Hamiltonian (jH) decouple 
in the freezing limit a — > oo, by Eq. ([7]) its eigenvalues 
are approximately given by 



a > 1 , 



(16) 



where Ef^ and £j are two arbitrary eigenvalues of Hgc 
and 7i, respectively. The asymptotic relation (|16p imme- 
diately yields the following exact formula for the partition 
function Z of the DAf-type PF spin chain ([9]): 



Z{T) = hm 



Z{aT) 



oo Zsc{aT) ' 



(17) 



where Z and Zsc are the partition functions of H and Hgc, 
respectively. Inserting the expressions for the partition 
functions Z and Zsc in the latter equation we shall obtain 
an explicit formula for the partition function Z of the 
chain 

In order to determine the spectra of the corresponding 
Hamiltonians H and Hsc in Eqs. (jH) and ([5]), following 



E 



a 



— {a-K,,) + ^-^{a-K,,] 



(18) 



where Kij and Ki are coordinate permutation and sign 
reversing operators, defined by 

{Kij f ){xi, . . . ,Xi, . . . ,Xj, . . . , xn) 

— / (^1 7 ■ ■ ■ ; -^j : • ■ ■ T ^i: • ■ ■ 7 "^N^ ; 

{Kif){xi, ...,Xi,.. .,xn) = f{xi, . . .,-Xi, . . .,xn) , 
and Kij ~ KiKjKij. We then have the obvious relations 



H = H'\ 

Hue = H'\ 



Kii-*l,Ki 



(19a) 
(19b) 



where e can take both values ±1. On the other hand, the 
spectrum of H' is easily computed by noting that this 
operator can be written in terms of the rational Dunkl 
operators oi type [13] 



Ji =9a.. +aXl 



X. 



-(l-if,,) + ^(l-xy) 



i = 1, . . . , N, as follows [3 

i 

where 

Eo = Na(^a{N 



^Xidxi + Ec 



1) 



(20) 



(21) 



(22) 



Since the Dunkl operators (|20p map any monomial Y[i x^' 
into a polynomial of total degree ni -|- • • • -f riAf — 1, by 
Eq. (j2ip the operator H' is represented by an upper tri- 
angular matrix in the (non-orthonormal) basis with ele- 
ments 

K^p\{xT, n= (ni,...,n^) e (NUlO})"^, (23) 



ordered according to the total degree |n| = ni -I- • • • -I- rijv 
of the monomial part. More precisely, 

(24) 

|m|<|n| 

where 

E'^ = a\n\+EQ (25) 
and the coefficients Cmn are real constants. 
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We shall now construct a basis of the Hilbert space of 
the Hamiltonian H in which this operator is also repre- 
sented by an upper triangular matrix. To this end, let us 
denote by S « (C™)®^ the Hilbert space of the su(m) 
internal degrees of freedom, and let 

|s) = |si,...,sjv), s, = -M,-M+l,...,M=^, 

be an arbitrary element of the canonical (orthonormal) 
basis in this space. Due to the impenetrable nature of 
the singularities of the Hamiltonian (j4|) , its Hilbert space 
is the set L/^{C) (g) S of spin wave functions square inte- 
grable on the open set C which vanish sufficiently fast on 
the singular hyperplanes Xi±Xj=Q,l^i<j^N. It 
can be shown, however, that H is equivalent to its nat- 
ural extension to the subspace of Lq(R^) (g) S consisting 
of spin wave functions antisymmetric under particle per- 
mutations and symmetric under sign reversals of an even 
number of coordinates and spins. (This is essentially due 
to the fact that any point in not lying on the singular 
subset Xi±Xj=0, l^i<i^N, can be mapped in 
a unique way to a point in C via a suitable element of 
the Djv Weyl group, which is generated by coordinate 
permutations and sign reversals of an even number of co- 
ordinates [1^.) We can therefore assume without loss of 
generality that the Hilbert space of H is the closure of 
the subspace spanned by the functions 

V'^,,(x)=A^(^„(x)|s)), e = ±l, (26) 

where A"^ denotes the projector on states antisymmet- 
ric under simultaneous permutations of spatial and spin 
coordinates, and with parity e under sign reversals of 
coordinates and spins. The latter functions are linearly 
independent, and hence form a (non-orthonormal) basis 
of the Hilbert space of -ff, provided that the quantum 
numbers n and s satisfy the following conditions: 

(i) ni'^ ■ ■ ■ riM- 

(ii) Si > Sj whenever = rij and i < j. 

(iii) Sj ;^ for aU i, and > if (-1)"' = 

The first two conditions are a consequence of the anti- 
symmetry of the states (j26p under particle permutations, 
while the last condition is due to the fact that V'n s miist 
have parity e under sign reversals. It should be noted 
that the Hilbert space V of the Hamiltonian H just de- 
fined can be written as the direct sum 

V = V+®V_, (27) 

where the subspace is the closure of the span of the 
basis vectors ■0n s(^)- Within each subspace V^, a partial 
ordering among these basis vectors may again be defined 
by the degree |n|. We shall now show that the Hamil- 
tonian H is represented by an upper triangular matrix 
in this basis (and thus by a direct sum of two upper tri- 
angular matrices in the total Hilbert space V). Indeed, 
since KijA'' = -SijA" and KiA'' = eSiA", it follows that 
HA" = H'A". Using this identity, Eq. and the fact 



that H' obviously commutes with A*^, we have 

^V'^,s=^V^,s = A^(Wn)|s)) 

= K^n,s+ E Cmn^m,s- (28) 
m| < I n 

Suppose now that both n and s satisfy conditions (i)- 
(iii) above, so that V'n,s belongs to the basis of under 
consideration. Although a given pair of quantum num- 
bers (m, s) in the r.h.s. of the previous equation need not 
satisfy these conditions, it is easy to see that there is a 
permutation TTm such that m' = 7rm(ni) and s' = 7rm(s) 
do satisfy (i)-(iii). Since V'm s differs from the basis vec- 
tor 'ip^, g, at most by a sign, and |m'| — |m| < |n|, our 
claim follows directly from Eq. Moreover, the latter 
equation and Eq. ((25)) imply that the eigenvalues of the 
spin Calogero Hamiltonian ^ are given by 

=a|n|+^o, (29) 

where e — ±1 and n, s satisfy conditions (i)-(iii) above. 
Since the numerical value of ^ is independent of s 
and e, the energy associated with a quantum number 
n will be highly degenerate in general. For any given 
n, this degeneracy factor d„ can be found by counting 
the numbers of independent spin states |s) satisfying 
conditions (ii) and (iii) for each case e = -1-1 and e — 
— 1, and finally taking the sum of these two numbers. 
Explicit expressions for such degeneracy factors will be 
given shortly when computing the partition function of 
the model. 

It is important at this point to elucidate the connection 
between the Hilbert spaces of the D^r-type spin Calogero 
model and its BCn counterpart. The key fact in this re- 
spect is that the Dpf Hamiltonian ^ does not depend 
on the discrete parameter e. Consequently, as shown in 
Eq. (pS)) . we can use both projectors A+ and A~ for con- 
structing the Hilbert space. On the other hand, since e 
appears explicitly in the Hamiltonian of the BCn spin 
Calogero model ([T]), for any given value of e only the 
corresponding projector A*^ can be used to construct the 
Hilbert space [2g|. Moreover, when 6 = this Hilbert 
space is essentially the subspace Vg of V in Eq. ([77)1 . 
Thus the presence of e in the Hamiltonian of the BCn 
spin Calogero model effectively introduces a "chirality" 
in this system. By Eq. (|77|) . the Hilbert space of the 
Dn spin Calogero model is simply the direct sum of the 
two Hilbert spaces associated with two BCn models with 
opposite chiralities (and 6 = 0). 

Turning next to the scalar Hamiltonian Hgc, in view 
of Eq. (jl9bp we now need to consider scalar functions of 
the form 

V'J.(x)=A^(/.„(x), e = ±l, (30) 

where A| is the projector onto states symmetric with re- 
spect to permutations and with parity e under sign re- 
versals. In fact, we can take as the Hilbert space of Hsc 
the space of symmetric functions in Lq(M") with even 
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parity with respect to an even number of coordinate sign 
reversals. In other words, the Hilbcrt space of i?sc is the 
direct sum of its two subspaces = AgLQ(E"), whose 
elements have parity e under sign reversals. The func- 
tions (|30p form a (non-orthonormal) basis of the corre- 
sponding subspace provided that either rii = 2ki for 
all i (for e = 1), or rii = 2ki + l for all i (for e = —1), with 
ki ^ • • • ^ fcjv in both cases. Just as before, if for each 
e = ±1 we order the basis functions ^/'n(x) according to 
the degree |n|, the matrix of the scalar Hamiltonian 7?sc 
in the basis ([50]) is expressed as a direct sum of two up- 
per triangular matrices, with diagonal elements also 
given by the r.h.s. of ((29)) . However, due to the absence 
in this case of spin degrees of freedom, the degeneracy 
factor of every quantum number n is one. Note also 
that from Eq. (|^^ . its analogue for the energies of the 
scalar Hamiltonian, and the freezing trick relation ([16]), 
it follows that all the energies of the spin chain ^ are 
integers. 

Let us next compute the partition functions Zgc and Z 
of the models ([5]) and ^ . To begin with, from now on we 
shall drop the common ground state energy Eq in both 
models, since by Eq. (fT7| it does not contribute to the 
partition function Z. With this convention, the partition 
function of the scalar Hamiltonian Hsc is given by 

Z,,(aT) = (l + g^) 9^'"' 

where q — e~^/('^^^). The latter sum is easily recognized 
as the partition function 

of the scalar Calogero model of BCn type evaluated in 
Ref. ;2J]. We thus have 

Z,,(aT) = (l + g^)Zif)(aT) 

=(i-g^)-in(i-'?'^)''- (31) 

We are now ready to compute the partition function of 
the spin Hamiltonian iJ in Eq. As for the BCj^ 

model [1^, it is convenient to deal separately with the 
cases of even and odd m. 



A. Even m 

When m is even, condition (iii) above simplifies to 
(iii') Si > for aU i. 

By Eq. ([29)) . after dropping Eq the partition function of 
the Hamiltonian ^ can be written as 

Z{aT)= J2 ^^"9'"'' (32) 

ni ^ ■ ■ -^n^ ^0 



where dn is the spin degeneracy factor associated with 
the quantum number n. Writing 

n ^ (^ki, . . . ,ki, . . . ,kr, . . . ,kr), /ci > • • • > fc^ > 0, 

(33) 

and using the conditions (ii) and (iii'), we have 

dn^2f[M^\^2d{,.), V^{u^,...,Vr), (34) 

where c?(f) is the corresponding degeneracy factor for the 
BCm type of spin Calogero model ([1]) with even m, and 
the factor of 2 is due to the two values taken by e in 
Eq. Note that YJi=i^i = so that the multi- 

index u can be regarded as an element of the set Vn of 
partitions of N (taking order into account). With the 
previous notation, Eq. p2p becomes 

Z{aT)^2Y, d{u) J2 9'^^""'" 

= 2Z(B)(aT), (35) 

where 

r TV- 

u^Vn j = l i=l 

is the partition function of the su(m) spin Calogero model 
of BCn type with even m, cf. [2^. From Eqs. (fT7)) . 
(1^ . (|55|) and the latter expression we finally obtain the 
following explicit formula for the partition function of the 
su(m) PF chain of _Djv type in the case of even m: 

N-l jv- 

z(r)=2n(i-a E ^('^) n (36) 

i=l v^Vn j=l 

where i^v) = r is the number of components of the multi- 
index u. The latter equation can be also written as 

N-i '<|r]Y N-e{i.) 

Z{T) = 2Y[il + q^) J2 d{v)q^-^ ' n ^-^'')^ 

(37) 

where the positive integers N'^ are defined by 

{N[,...,N'^_,^,)] = {!,..., iV-l}-{7Vi,...,iV,(,)_i}. 

Note also that from the freezing trick relation (fT7|) . its 
analogous for the BCn models, and Eqs. ([3T|) - ([35|) one 
easily obtains the identity 

Z(T) = 2(l + g^)-iz(B)(T) (evenm), (38) 

where Z^^\T) is the partition function of the su(m) PF 
chain ^ of BCn type. 
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For the simplest case of spin 1/2 chain, we have Vi = 1 
for aU z, and therefore l{i>) = iV, d{u) — 1 and Nj = j, 
so that Eq. ([37|l simpHfies to 



AT-l 



2(T) = 2(7^^(^-1) [| (1 + , m = 2. (39) 

1=1 

Thus, for spin 1/2 the spectrum of the chain ([9]) is given 
by 

£j = ^NiN-l)+j, J =0,l,...,i7V(7V-l), (40) 

and the degeneracy of the energy £j is twice the number 
Qn-iU) of partitions of the integer j into distinct parts 
no larger than iV — 1 (with (3jv_i(0) = 1). 



where Zjf denote the partition functions of the su(m) 
spin Calogero models of BCn type H]) with odd m and 

e = ±1. Using the expressions of derived in Ref. [Ill 
we finally obtain 



-{N+N,) TT 9 



nr 



(44) 



j=s+l 



Substituting the previous expression and (PT]) into (|17p . 
we immediately deduce the following explicit formula for 
the partition function of the su(r7i) PF chain of Dn type 
for odd m: 



B. Odd m 

Let us consider now the case of odd m. As for the 
BC'n chain, in this case it is convenient to slightly modify 
condition (i) above by first grouping the components of 
n with the same parity and then ordering separately the 
even and odd components. In other words, we shall write 
n = (oe, Ho), where 



N-l 



z(r) = (i-g^) nii-?'^) E E^«('^)« 



nr 



n,i 



2Ni 



(45) 



Equivalently (cf. Eqs. ([3T|) and (|43l) ) 



no = (2fci, . . . ,2fci, . . 



,2fc.). 



no = ( 2fc,+i + 1, . . . , 2fc,+i + 1, . . . , 2A:^ + 1, . . . , 2fc^ + 1 ) , 
and 

ki > ■ ■ ■ > ^ 0, /cs+i > • • • > fcr ^ . 

The spin degeneracy factor is now 

dn = d-iu) + dt{u)=ds{u), (41) 

where df{v>) is the number of independent spin states |s) 
satisfying conditions (ii) and (iii) with e = ±1, namely 
(cf. [2i, Eq. (28)]) 



2 



n 



7n~€ 
2 



(42) 



Calhng 



Nj = ^ i^i, j = s + 1, . . . ,r , 



Z{T) = (1 + g^)-i(zf ^(T) + Zi^^(r)) (odd to) , 

(46) 

where (T) are the partition functions of the su(m) 
PF chains ([3]) of BCn type for odd to. Note that the 
latter formula is also valid for even to, since in that case 
= Z^_^^ = Z^^\ In fact, Eq. dH]) can be used to 
verify that the expression (|45p for the partition function 
of the su(to) PF spin chain oi type is a polynomial 
in q, as should be the case for a finite system with inte- 
ger energies. To this end, recall from Ref. [s^l that the 
partition function zi^'' can be written as 



N 



K=0 
N 

X J] + 

i=K+l 



r(A) 
■'N-K 



{q;^), (47) 



where Z^^j^{q; ^^^^^) is the partition function of the 
su(2^^^^) PF spin chain of An type with N — K particles, 
and 



and proceeding as before, we obtain 



E 2iy.fe, E 1^.(2^.-1-1) 



,= 1 + ! 



'(ar) + zi^^(aT), 



(43) 



N 
K 



N 



{q)K{q)N-K 



{q),^Y[{l-q^) 



It can be shown that both the g-binomial coefficient [ ^ ] 



and the partition function Z^^j^ are polynomials in q, 
cf. Refs. [S^, Since all the terms in the sum in the 
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r.h.s. of Eq. p7|) contain a factor of 1 + g except for 
K — N , the partition function Zc ' can be expressed as 

(T) = (1 + q'')VM + g^^^-^)?* (^"^^ , 



where 



Af-l 



^=0 



r(A) 



is a polynomial in g. Inserting the latter equations 
into (pSj) we immediately conclude that 



N 



N-1 



K=0 i=K 

is a polynomial in g, as claimed. 



N' 
K 



7(A) j m-l \ 
■'N~K\'i' 2 



(48) 



IV. STATISTICAL ANALYSIS OF THE 
SPECTRUM 

In this subsection we shall take advantage of the ex- 
plicit expressions for the partition function of the su(77i) 
PF chain of Dn type ^ just derived to check that its 
spectrum shares the global properties of those of other 
spin chains of Haldane-Shastry type mentioned in the 
Introduction. In practice, in order to compute the spec- 
trum for given values of N and m it is more efficient to 
use Eq. (|48l) for odd m and its analog for even m 



N-l 



Z(T) = 2z(f)(9;f) 



(49) 



j=i 



obtained from Eq. §^ using Eq. (31) in Ref. to- 
gether with the explicit expression 



E 



. E M,{M,-1) 



(9) 



K 



A/iH hM„ 



Ml 



(9) 



derived in Ref. [30|. With the help of the previous for- 
mulas it is possible to determine the chain's spectrum 
for relatively large values of N and m; for instance, using 
Mathematica"""^ on a personal computer it takes less 
than 10 seconds to evaluate the partition function in the 
case iV = 50 and m = 3. 

In the first place, our calculations of the spectrum for 
a wide range of values of m and N show that the ener- 
gies of the Dff chain ([9]) form a set of consecutive inte- 
gers, as is the case for all the previously studied (non- 
supersymmetric) rational chains, of both An and BCn 
type [H, HI] . As to the (normalized) level density 



fin 



di 5{£ — £i 



(50) 



where £1 < ■ ■ ■ < £l are the distinct energy levels and 
di is the degeneracy of £i, we have verified that when 
N is sufficiently large it can be approximated with great 
accuracy by the Gaussian law 



5(f) = 



1 



27r(T 



(51) 



with parameters fi and a given by the mean and standard 
deviation of the chain's spectrum. Since the energy levels 
are consecutive integers, this means that 



g{£,) {N » 1) . 



(52) 



As an illustration, in Fig. [T] we have plotted both sides 
of the latter equation in the case m = 2 and N = 20. 



0.015 



0.010 



0.005 




200 



250 



300 



350 



FIG. 1: Plot of the Gaussian distribution (|5ip (continuous 
red line) versus the l.h.s. of Eq. (|52p (blue dots) in the case 
m — 2 and A'^ = 20. The root mean square error (normalized 
to the mean) of the adjustment is 3.01 x 10^^. 



In view of the approximate relation ([5^ . it is of inter- 
est to evaluate the mean and standard deviation of the 
energy in closed form for arbitrary values of N and m. 
This can be done in essentially the same way as for the 
BCm chain using the formulas for the traces of the 
spin operators Sij, Si and Sij in Ref. [2^. Indeed, setting 



the mean energy is given by 



H = trH = (1 + — ) '^{h,j + hij) 



The sum in the r.h.s. of the previous equation is clearly 
half the maximum energy £niax of the Hamiltonian 
so that by Eq. ((6T|) we have 



2 V TO/ 



(53) 
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Similarly, the variance of the energy is given by 



One then defines the normalized spacings 



where y is the parity of to, and we have used Eq. (A6) in 
Ref. [il. From Eqs. (A8), (A9) and (A12) of the latter 
reference with /3 = 0, one easily obtains 



_ 1 / 1 
36 V m? 



N{N~l){AN+l)-- 



1 



■{\~p)N{N-l). 

(54) 

With the help of the above expressions for and cr, 
we can show that Eq. is not incompatible with the 
fact that the level densities of the three chains Ti. and 
Ti^^) with e = ±1 are approximately Gaussian for large 
N. Indeed, writing (|^ as 



ll + q''')Z{T)^zf\T) + Zf>{T) 



r(B) 



(55) 



we see that the l.h.s. of ([55)1 represents the superposition 
of the spectrum of the Dm chain ([9]) and its translation 
by TV, whose level density tends to the sum of the Gaus- 
sian g{S) in ([5T|) and its translate g{£ — N) as N oo. 
But in this hmit we have iV < tr = 0{N^^'^), so that 
g[E) + g{E: + TV) ~ 2g{£). Similarly, the r.h.s. of Eq. ^ 
is the partition function of the superposition of the spec- 
tra of the chain Hamiltonians ([3]) with e = ±1, whose 
level density for large N is approximately the sum of two 
Gaussians with the same standard deviation as (|5ip and 
mean equal to -I- ^ (l — , cf. Ref. [28i|. Hence as 
— > oo the level density of the r.h.s. of ([55]) is approxi- 
mately given by 



»(^-f('-^))-K^-f('-^))-^»<^). 

as the l.h.s. 

Let us consider now the distribution of the spacings 
between consecutive levels in the "unfolded" spectrum. 
Recall [ill, to begin with, that the unfolding of the levels 
£i of a spectrum is the mapping £i i— > 77,; = tjiSi), where 
r]{£) is the continuous part of the cumulative level density 



-N 



The unfolding mapping makes it possible to compare dif- 
ferent spectra in a coherent way, since the unfolded spec- 
trum {f^i}^! can be shown to be uniformly distributed 
regardless of the initial level density. In our case, by 
the above discussion we can take ?/(5) as the cumulative 
Gaussian density (|5ip . namely 



1 ^ 



1 + erf 



( V2<J 



(56) 



im+i - Vi)!^-, 



i = 1, . . . ,L - 1 , 



where A = {rji^ — r]i)/{L — 1) is the mean spacing of 
the unfolded energies, so that {si}fs^^ has unit mean. 
According to a well-known conjecture of Berry and Ta- 
bor, for a quantum integrable system the density p{s) 
of normalized spacings should be given by Poisson's law 
p{s) = c^**. By contrast, for a system whose classi- 
cal counterpart is chaotic, it is generally believed that 
the spacings distribution follows instead Wigner's law 
p{s) — {tts/2) exp(— 7rs^/4), typical of the Gaussian en- 
sembles in random matrix theory [JH . 

We shall now see that the spacings distribution of the 
PF chain of Dn type © follows neither Poisson's nor 
Wigner's law, as is the case for all spin chains of HS 
type studied so far pT?, [Hi IH, Es] . More precisely, 
we will show that the cumulative spacings distribution 
P{s) = Jp p(s')ds' is approximately given by 



P{.s) 1 - 



log 



(57) 



where Smax is the maximum spacing. In fact, as proved 
in Ref. [2^ , the previous approximation necessarily holds 
for any spectrum £min = £1 < ■ ■ ■ < £l ^ £miix satisfying 
the following conditions: 

(i) The energies are equally spaced, i.e., £i+i — £i — 5£ 
for i = 1 , . . . , L — 1 . 

(ii) The level density (normalized to unity) is approxi- 
mately given by the Gaussian law (|5ip . 



(iii) £„ 



fi, fj,-£„ 



> cr. 



(iv) £niin and £max are approximately symmetric with 
respect to namely l^min -I- fmax - 2/i| < fmax - 

Moreover, when these conditions are satisfied the maxi- 
mum spacing can be estimated with great accuracy as 

£ — £ ■ 



(58) 



It should also be noted that Eq. ([57]) is valid only for 
spacings s € [so, Smax], where 



So 



(59) 



is the unique zero of the r.h.s. of (j57p (the inequality 
in (|59|) follows easily from condition (iii) and Eq. ((58| V 
We shall next check that conditions (i)-(iv) above are 
indeed satisfied by the spectrum of the chain ([9|) when 
A^ ^ 1. In fact, we already known that conditions (i) 
(with 8£ = 1) and (ii) hold. In order to verify condition 
(iii), we first need to compute the maximum and mini- 
mum energies fmax and fmin- The maximum energy is 
clearly 



^max — 2 ^ ^ |^(^j ^j) 



+ (6+ej)" 



(60) 
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whose corresponding eigenvectors are the spin states 
symmetric under permutations and with parity ±1 under 
spin reversals. Since fmax is independent of m, it is most 
easily computed for the spin 1/2 chain, whose spectrum 
is explicitly given in Eq. (|40)l . We thus obtain 



fmax = N{N - 1) . (61) 

As to the minimum energy, Eq. (j55[) implies that 
£ ■ -minff'^' f(^) ) 

<--min — mm V'-'min,- ' '-'min, + i ' 

where the minimum energies f ^ of the BC n chain ([3]) 



were computed in Ref. [28l. From Eqs. (B1)-(B2) of the 

'mm,-' 



latter reference it easily follows that f/^J^j^ ^ E^^^ 
that 



N 



N 



^min — ^ I 1 



m 



m/ Zm 

X {l-mp9{2l-m-l)), (62) 



with 



/ = N mod — (1 
2 ^ 



P) 



From Eqs. ^ and ^ it immediately follows 

that (fmax — m)/"" and (^min — A*)/^ both 0{N^^'^) as 
iV — > oo, so that condition (iii) is also satisfied. Finally, 
from the latter equations it also follows that £min+'S^max — 
2fi is at most 0{N) while £,nax — £min = 0{N'^), which 
proves condition (iv). 

The previous argument shows that the cumulative 
spacings distribution of the chain ^ should be well 
approximated by the r.h.s. of Eq. ([57)1 when N is suf- 
ficiently large. We have verified that ([57]) is indeed in 
excellent agreement with the numerical data for many 
different values of N and m. For instance, in the case 
TV = 20 and to = 2 presented in Fig. [5] the root mean 
square error (normalized to the mean) of the adjustment 
of P(s) to the r.h.s. of Eq. §7} is 1.03 x 10"^^ and this er- 
ror decreases to 4.69 x 10^* when N = 100. It should be 
stressed that the approximation (j57p contains no free pa- 
rameters, since the maximum spacing Smax is completely 
determined as a function of N and to by Eqs. 
(pT|) and (p^ . In fact, from the latter equations it im- 
mediately follows that for large N the maximum spacing 
is asymptotically given by 



y m + I 



(63) 



as for the (non-supersymmetric) PF chains of BCn 
type M- 




0.5 1.0 1.5 2.0 2.5 3.0 



FIG. 2: Cumulative spacings distribution P{s) and its ap- 
proximation (|57[) (continuous red line) for A'' = 20 and m = 2. 
For convenience, we have also represented Poisson's (green, 
long dashes) and Wigner's (green, short dashes) cumulative 
distributions. 



and its corresponding spin model 

,2 



can be studied in much the same way as their antifcrro- 
magnetic versions (HJ-dHl). Since now 



Hv = H'\ 



(65) 



we must replace the operator A' in Eq. ([^5]) by the pro- 
jector Ag onto states symmetric under simultaneous per- 
mutations of the particles' spatial and spin coordinates, 
and with parity e under sign reversal of coordinates and 
spin. Hence condition (ii) above for the new basis states 



Vij.,,^A^(0„(x)|s)), 



±1, 



should now read 

(ii') Si ^ Sj whenever rii = rij and i < j. 

As a result, the degeneracy factors d{v>) and ds{i^) in 
Eqs. and (|¥T|) should be replaced by their "bosonic" 
versions 



and 



THE FERROMAGNETIC CASE 



where 



The ferromagnetic spin chain of Dn type with Hamil- 
tonian 



^^3 



1 - 



1 - 



(64) 



1=1 



2 



n 

i=s+l 



ra — e 
2 



Therefore the partition function of the ferromagnetic 
su(m) PF chain oi type ([M]) is still given by Eq. 
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(for even to) or (l45l) (for odd to), but with d(f ) and ds{i') 
replaced respectively by dF{u) and dp^sii^)- 

On the other hand, the chains ([9]) and ([64|) are obvi- 
ously related by 

T^F + 7^ = 2 ^ [(6 - + (C. + 0)-'] = ^(A^ - 1) , 

(66) 

where we have used Eqs. (|60 |) -(|6T |) . Thus the partition 
functions 2 and Zp of Ti. and Tip satisfy the remarkable 
identity 



Zp(g)=g^(^-i)Z(g-i) 



(67) 



This is a manifestation of the boson-fermion duality 
discussed in detail in Refs. for the su(TO|n) supersym- 
metric HS spin chain, since the ferromagnetic (resp. an- 
tiferromagnetic) chain can be regarded as purely bosonic 
(resp. fermionic). For instance, using the latter identity 
and Eq. ([39]) we easily obtain the following expression 
for the partition function of the ferromagnetic spin 1/2 
chain: 



N-l 



Zf(T) = 2 [](l + g'), m = 2. 



(68) 



With the help of the duality relation ((67)) and the ele- 
mentary g-number identity 



K 



it is straightforward to derive the analogs of Eqs. (|48p and 
(gni) for the ferromagnetic chain CaUing Z^p{q; n) 



the partition function of the su(n) ferromagnetic PF 
chain of type A for K spins, given by [sO] 



K 



we obtain in this way 

JV-l 



Z(T) = 2Z^^)(g;^) \{{l + q^) 
for even to, and 



N N-l 
K=0 i=K 



^W-if.FWj 2 ) 



for odd TO. Finally, from the duality relation (I66p it 
clearly follows that the statistical properties of the spec- 
trum of Tip are identical to those of TL^ namely when N 
is large enough the level density is approximately Gaus- 
sian, and the spacings distribution follows Eq. l[57|) with 
great accuracy. 

Acknowledgments 

This work was partially supported by Spain's DGI 
under grant no. FIS2005-00752, and by the Com- 
plutense University and Madrid's DGUI under grant 
no. GR74/07-910556. 



[1] F. Calogero, J. Math. Phys. 12, 419 (1971). 
[2] B. Sutherland, Phys. Rev. A 4, 2019 (1971). 
[3] B. Sutherland, Phys. Rev. A 5, 1372 (1972). 
[4] M. A. Olshanetsky and A. M. Perelomov, Phys. Rep. 94, 
313 (1983). 

[5] F. D. M. Haldane, Phys. Rev. Lett. 60, 635 (1988). 

[6] B. S. Shastry, Phys. Rev. Lett. 60, 639 (1988). 

[7] A. P. Polychronakos, Phys. Rev. Lett. 70, 2329 (1993). 

[8] Z. N. C. Ha, Quantum Many-body Systems tn one Dimen- 
sion, vol. 12 of Advances in Statistical Mechanics (World 
Scientific, Singapore, 1996). 

[9] M. V. N. Murthy and R. Shankar, Phys. Rev. Lett. 73, 
3331 (1994). 

[10] A. P. Polychronakos, J. Phys. A 39, 12793 (2006). 
[11] H. Azuma and S. Iso, Phys. Lett. B 331, 107 (1994). 
[12] C. W. J. Beenakker and B. Rajaei, Phys. Rev. B 49, 7499 
(1994). 

[13] M. Caselle, Phys. Rev. Lett. 74, 2776 (1995). 

[14] N. Taniguchi, B. S. Shastry, and B. L. Ahshuler, Phys. 

Rev. Lett. 75, 3724 (1995). 
[15] P. J. Forrester, Nucl. Phys. B416, 377 (1994). 
[16] H. Ujino and M. Wadati, J. Phys. Soc. Jpn. 66, 345 

(1997). 



[17] T. H. Baker and P. J. Forrester, Nucl. Phys. B492, 682 

(1997) . 

[18] D. Bernard, M. Gaudin, F. D. M. Haldane, and 

V. Pasquier, J. Phys. A 26, 5219 (1993). 
[19] K. Hikami, Nucl. Phys. B441, 530 (1995). 
[20] B. Basu-Mallick and A. Kundu, Nucl. Phys. B509, 705 

(1998) . 

[21] Z. N. C. Ha and F. D. M. Haldane, Phys. Rev. B 46, 
9359 (1992). 

[22] K. Hikami and M. Wadati, J. Phys. Soc. Jpn. 62, 469 
(1993). 

[23] J. A. Minahan and A. P. Polychronakos, Phys. Lett. B 

302, 265 (1993). 
[24] B. Sutherland and B. S. Shastry, Phys. Rev. Lett. 71, 5 

(1993). 

[25] H. Frahm, J. Phys. A 26, L473 (1993). 

[26] A. P. Polychronakos, Nucl. Phys. B419, 553 (1994). 

[27] F. Finkel and A. Gonzalez-Lopez, Phys. Rev. B 72, 

174411(6) (2005). 
[28] J. C. Barba, F. Finkel, A. Gonzalez-Lopez, and M. A. 

Rodriguez, Phys. Rev. B 77, 214422(10) (2008). 
[29] A. Enciso, F. Finkel, A. Gonzalez-Lopez, and M. A. 

Rodriguez, Nucl. Phys. B707, 553 (2005). 



12 



[30] B. Basu-Mallick, H. Ujino, and M. Wadati, J. Phys. Soc. 

Jpn. 68, 3219 (1999). 
[31] B. Basu-Mallick and N. Bondyopadhaya, Nucl. Phys. 

B757, 280 (2006). 
[32] J. C. Barba, F. Finkel, A. Gonzalez-Lopez, and 

M. A. Rodrfguez (2008), Nucl. Phys. B, in press 

(doi: 10. 1016/j .nuclphysb.2008.08.014) . 
[33] J. C. Barba, F. Finkel, A. Gonzalez-Lopez, and M. A. 

Rodriguez, Europhys. Lett. 83, 27005(6) (2008). 
[34] M. V. Berry and M. Tabor, Proc. R. Soc. Lond. A 356, 

375 (1977). 

[35] E. Corrigan and R. Sasaki, J. Phys. A 35, 7017 (2002). 
[36] S. Ahmed, Lett. Nuovo Gimento 22, 371 (1978). 



[37] C. F. Dunkl, Commun. Math. Phys. 197, 451 (1998). 
[38] F. Finkel, D. Gomez-Ullate, A. Gonzalez-Lopez, M. A. 

Rodriguez, and R. Zhdanov, Nucl. Phys. B613, 472 

(2001). 

[39] J. E. Humphreys, Introduction to Lie Algebras and Rep- 
resentation Theory, Graduate Texts in Mathematics 9 
(Springer- Verlag, New York, 1972). 

[40] J. Ciglor, Monatsh. Math. 88, 87 (1979). 

[41] T. Guhr, A. Miillcr-Grochng, and H. A. Wcidenmiiller, 
Phys. Rep. 299, 189 (1998). 

[42] B. Basu-Mallick, N. Bondyopadhaya, K. Hikami, and 
D. Sen, Nucl. Phys. B782, 276 (2007). 



